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Abstract. - We use the weak coupling renormalization group method to examine the interplay 
between charge-density- wave and s-wave superconducting orders in a quasi-one-dimensional model 
of electrons interacting with acoustic phonons. The relative stability of both types of order is 
mapped out at arbitrary nesting deviations and Debye phonon frequency ujd- We singled out 



a power law increase of the superconducting Tc 



from a quantum critical point of charge- 



density-wave order triggered by nesting alterations. The results capture the key features shown by 
the proximity between the two types of ordering in the phase diagram of the recently discovered 
Perylene based organic superconductor under pressure. The impact of Coulomb interaction on 
the relative stability of the competing phases is examined and discussed in connection with the 
occurrence of s-wave superconductivity in low dimensional charge- density- wave materials. 
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Introduction. — The recent observation of supereon- 
uctivity (SC) in the Perylene based organic conductor 
Per2[Au(nint)2] [I], brings once again into focus the possi- 
ble role played by a charge density- wave (CDW) instability 
ih the mechanism of onset of SC in quasi-one-dimensional 
(quasi-lD) electron systems. This worh makes use of the 
Renormalization group (RG) technique to analyze the in- 
terplay between these two phases in systems of electrons 
doupled to lattice vibrations. Besides their relevance for 
materials showing a CDW-SC proximity, the results also 
dddress the issue of quantum criticality associated with in- 
terfering orders for models of electrons coupled to bosonic 
^citations in low dimensions 

' Per2[Au(mnt)2] is a member of the two-chain charge 
transfer salts series Per2[M(mnt)2], where M= Pt, Pd, 
Au, .... These organic salts are made of Perylene and 
Dithiolate flat molecular complexes that pile up as seg- 
regated stacks well described by a quasi-lD electronic 
structure @]. In normal pressure conditions, the Pery- 
lene chains undergo a metal-insulator transition due to 
the formation of a Peierls lattice distorted state driven by 
a CDW superstructure [5]. For the Per2[Au(mnt)2] com- 
pound, only the Perylene stacks are electronically active in 
the CDW transition, which takes place at Tcdw — 12 K 
at ambient pressure [HJ. This is a relatively low temper- 
ature scale likely to be vulnerable to nesting alterations 



of the Fermi surface by the application of pressure. This 
is supported by the suppression of the insulating state 
under 5 kbar of pressure, which turns out to be also crit- 
ical to the onset of SC at Tc ~ 300 mK [T], hinting at a 
direct part played by CDW correlations in the enhance- 
ment of Cooper pairing. The sequence of states thus ob- 
tained is reminiscent of the competition between CDW 
and SC orders found in some quasi-lD transition-metal 
trichalcogenides materials [7j. The pattern is also akin to 
the quasi-lD Bcchgaard salts series [(TMTSF)2X], where 
a spin-density-wave (SDW) state is known to be followed 
by superconductivity under pressure [SUlOj. In the latter 
case the application of the RG method to electron mod- 
els with momentum dependent repulsive interactions has 
demonstrated how nesting deviations can control the in- 
terference between density-wave and non s-wave Cooper 
pairings and reproduce the sequence of phases displayed 
by the Bechgaard salts under pressure [TTH13] . 

However, at variance with SDW systems where the di- 
rect Coulomb term dominates the scene of interactions, 
the electron-phonon coupling plays an essential role in 
lattice distorted CDW systems. Electron-electron interac- 
tions induced by the exchange of acoustic phonons are dy- 
namically governed by the Debye energy scale lod, which 
is much smaller than the Fermi energy and often close to 
the energy scale of CDW order found in molecular conduc- 
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tors [5]. This introduces retardation in interactions, which 
besides interchain hopping and nesting aherations, mod- 
ifies in a non trivial way the interfering many-body pro- 
cesses that are Hnked to density-wave and Cooper pairings 
in every order of perturbation theory. This difficulty has 
been well established in the past, requiring to go beyond 
the habitual scheme of approximations such as mean-field 
and RPA-likc approaches that arc known to single out 
one pairing channel to the detriment of the other [T4HT6] . 
In the one-dimensional case, a weak coupling solution to 
this problem has been found in the framework of the RG 
method [T71[TH]. Recent progress along these lines has 
shown that this approach is well suited to simultaneously 
account for both pairing processes in the determination 
of ground states in electron-phonon systems at arbitrary 
phonon frequency. 

In this paper the RG method is extended to a electron- 
phonon model in the quasi-lD case and at finite temper- 
ature. The temperature scales Tcdw a-nd Tc for the in- 
stabilities of the metallic state against the formation of 
CDW and SC orders arc determined for arbitrary phonon 
frequency uid and nesting deviations parametrized by the 
next-to- nearest neighbor interchain hopping t'^ . The main 
results of the present work are outlined in the phase dia- 
gram of Fig. [TJ For small Tcdw weakens and under- 
goes a quantum-classical crossover as cj/^ is raised and goes 
beyond the adiabatic scale Tq for CDW ordering. When 
nesting distortion attains some threshold t'^ , Tcdw is crit- 
ically reduced and at non zero oj^i, an SC instability takes 
place in the s-wavc channel only. In the adiabatic limit, 
Tcdw defines a quantum critical point at t'^ , from which 
an anomalous power law increase of Tc with the phonon 
frequency takes place. Along realsitic pressure paths in 
the [t'^^LOr)) plane, the model phase diagram follows the 
leading features displayed by Per2[Au(mnt)2]. The impact 
of the repulsive Coulomb interaction on the structure of 
the phase diagram is explored alongside the predisposition 
of electron-phonon driven CDW systems to show s-wave 
superconductivity. 

The model and the renormalization group equa- 
tions. — We consider a non half-filled two-dimensional 
electron system consisting of N±_ chains of length L with 
the electron spectrum Ep{k) = vp{pk — kp) + e±{k±), 
where e_L(fcj^) = ~2t±cosk± — 2t^cos2A:^. Here p = ± 
refers to right and left moving electrons along the stacks, 
vpikp) to the parallel Fermi velocity (wave vector), and 
tj^ to the interchain hopping integral. In the quasi-lD 
case, we have t±_ <C Ep = vpkp, where Ep = i?o/2 
is the Fermi energy taken as half the band width. The 
next-to- nearest neighbor transverse hopping , which de- 
scribes nesting deviations, is kept small compared to t_\_. 
The following calculations are carried out for the typical 
values Ep = 15i_L = 3000 K. In the framework of the 
Su-Schrieffer-Heeger model [H], the electrons are linearly 
coupled to parallel acoustic phonons. These modes being 
harmonic, this is equivalent in the Matsubara formalism of 




Fig. 1: Phase diagram of the quasi- ID electron-phonon model 
in the {t'^,ujD) plane. Here To is the CDW ordering temper- 
ature of the adiabatic ujd and perfect nesting t'^ ^ G 
limit. 



the partition function to consider a frequency dependent 
electron-electron interaction. In the g-ology picture, the 
bare interaction reads 



1 + {uJni - U!ru)^/i^j 



(1) 



which for a non half-filled band splits as a backward (i = 1) 
and forward {i = 2) scattering amplitude between p ^ 
+ and — moving carriers; here the momentum-frequency 
variables k = (fc_L,w„) satisfy the conservation rule k^ + 
ki = ki + k2. In the following calculations we shall use the 
bare initial amplitudes (normalized by nvp) gf^ = —0.20 
for the backscattcring part, while = for the forward 
scattering at vanishing momentum transfer - the latter 
acquiring finite values following the RG transformation. 

The RG transformation of the coupling constants re- 
sults from the successive integration of electronic degrees 
of freedom in the outer energy shell ±Eo{£)d£/2 above 
and below the Fermi surface for all Matsubara frequencies. 
Here Eq{£) ~ Eqc"^ is the scaled bandwidth at step £ > 0. 
In the momentum-frequency RG scheme adopted here at 
finite temperature, each constant energy sheet from the 
Fermi surface is divided into 12 patches, each defining a 
particular fc^ in momentum space, while a discrete set of 
= 14 fermion Matsubara frequencies a;„ (—7 < n < 6) 
is retained along the frequency axis. At finite temperature 
this represents a good compromise between exacting com- 
puting time and reproducing the results known for either 
the non retarded case in quasi-one dimension |12j or the 
electron-phonon problem in one dimension |17| . 

At the one-loop level, the backward and forward scat- 
tering amplitudes obey the fiow equations 



dtgi{ki,k2,kz,ki) 



1 

2^ 



dkj_Ip{k±,qp) 



f-p{9i{ki,k, kp, ki)gi{kp, ^2,^3,^)) 
+ f-P,v{g2{ki,k, ki, kp)gi{kp, ^2,^3, k)) 
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and 



+ ep,v{gi{ki,k, kp, ki)g2{kp, h, k, h)) 
+ ^ y" dk^Icik_i_,qc) 
<^c{gi{ki,k2, k, kc)g2{k, fee, ^4, ^3)) 
+ ec(.92(fci, ^2, kc, k)gi(k, kch.h)) 



deg2{ki,k2,k3,k4) = ^ J dk±Ip{k±,q'p) 

X ep^i{g2[ki,k,ks„k'p)g2[k'p,k2,k,ki)) 
1 



(2) 



27r 



dk±Icik±,qc) 



ec (ffi (^1 ,k2,k, kc)gi {k,kc,ki,k3)} 

+ ec(52(fci,fc2,^c,fc).92(fc, ^0,^3,^4)) • (3) 



These consist of closed loop (ep = —2), vertex correc- 
tions (ep,t, = 1) and ladder {epj = 1) diagrams of the 
2fcp electron- hole (Peierls) pairing, which combine with 
the ladder diagrams (ec = —1) of the electron-electron 
(Cooper) pairing. Here kp = k+qp ,k'p = k+q'p and kc = 
—k + qc, where qp^c = {q±p,c,'-^p,c) corresponds to the 
Peierls qp = ki—k^, q'p = ki — k^ and Cooper qc = ^2 + ^1 
variables. In the above equations, each diagram contains 
a frequency convolution of the form gigj^c.p, which 
has been decoupled as {gigj)^^^ C.c.p at finite temper- 



ature. Here (• • • ) 



^„ • • ■ , stands as an average 



of the couplings over the internal loop frequency variable, 
whereas the (. derivative of the Cooper and Peierls loops 



p,c 



E 



-t-00 



^ £p,c is evaluated exactly to give 



Ip,c{k^,qp^c) = e[|i?oW/2 + i^Ap,c|-£;oW/2] 



1 

^4 



AT 

iEo{i) + iyAp,c)Eoii) 
{Eoi£) + iyAp^cr+oj%^c' 



4T 



(4) 



where Ap = -e{kj_)-e{k±+q±p), Ac = -e(/cj_)+e(fc^ + 
q±c), and Q{x) is the step function (8(0) = ^). 

A singularity can occur in either the Peierls or the 
Cooper scattering channel, which is indicative of an in- 
stability of the metallic state towards long-range or- 
der at the temperature T^. The nature of the or- 
dering state is provided by the singularity of the 
static and normalized response function TrvpXfj.i'i'^) = 
(27r)~-^// dk±di {z'^{k))Ip_c{k±, q^ ^) at the wave vector 
qO.^^ = (2fcp, tt) for = CDW and qgc = for fi = SC. 
The response vertex parts are governed by the equa- 
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Fig. 2: Typical low temperature dependence of the charge- 
density- wave (xcDw: squares) and s-wave superconducting 
(Xsc: triangles) normalized susceptibilities in the a) CDW and 
b) SC ordering sectors of Fig. |4l 



tions 



dezcBwik + qp) = — J dk^Ip{k^,qp)zcDwik + <?p) 

x{[epgi{k'+fp,k,k + q°p,k') 

+ ep,vg2{k'+ fp,k,k',k + 4)]), 



(5) 



for the /i =CDW response {qp ~ (tt, O)) and 



dizsci-'k + q'c) = dk'^Ic{k'^,fc)^sc{~~k' + fc) 
x{ec[gii-k'+q°c,k',~k + fc,k) 

+ g2{-~k' + (fc~k',k,~k + (fc)]). 

(6) 

for the static (s-wave) ^ =SC response (g^ = 0). For the 
whole range of parameters covered by the present model, 
the finite temperature singularities only occur for either 
the CDW or s-wave SC susceptibilities (Fig. [2|). 

Results. — Let us first consider the instability of the 
metallic state as one moves along the phonon frequency 
axis at fixed t'^ (Fig. [T]). At perfect nesting t'^ = 0, the 
adiabatic limit wp) — > is characterized by a singular- 
ity signaling the occurrence of a Peierls instability at the 
temperature denoted Tq (~ 20K for the parameters cho- 
sen here). In this limit, only close loops contribute to the 
fiow of Eqs. (I2]|3l) and ([5]), a limit equivalent to the molec- 
ular field analysis of the Peierls instability of the metallic 
state. By increasing wp, both the vertex and ladder di- 
agrams are progressively unlocked and begin to mix and 
interfere with closed loops. In the pertinent temperature 
range T <C t_L where all the instabilities take place, the 
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Fig. 3: (a) Normalized Tcdw (square) and s-wave Tc (trian- 
gle) vs the phonon frequency ratio lud/2To at different t'^/t±. 
(b) Normalized Tcdw (square) and s-wave Tc (triangle) vs the 
nesting deviation parameter t'^/t± for different frequency ra- 
tios ujd/2To. 



transverse electronic motion is coherent. As a function 
of energy, the interference is then maximum in the one- 
dimensional part of the flow where Eo{£)/2 > t±, whereas 
for Eq{£)/2 < t± interchain hopping begins to be coherent 
and the interference becomes non uniform in momentum 
space and generates a k± dependence of the coupling con- 
stants pniT^ . As a result, Tcdw diminishes with increas- 
ing ujd. However, when the frequency reaches the classical 
Peierls scale uj}j{t'j_ = 0) ~ 2To, the decrease is more rapid 
and Tcdw undergoes a crossover toward a non adiabatic 
CDW regime where all diagrams of both pairing channels 
contribute and ultimately level off the reduction of Tcdw 
(Fig. [31-a) - a crossover analogous to the one found in the 
purely ID case at tj^ = [TOIITT] . 

A finite amplitude of the anti-nesting term t'j_ modi- 
fies Ap in which reduces all the Peierls diagrams 
but leaves those of the Cooper channel unchanged. Thus 
Tcdw first gets smaller (Fig. [21-a), then drops rapidly as 
the phonon frequency extends across uj^, which is also 
smaller, to finally attain the non adiabatic limit (Fig.[3]-a). 
This reduction carries on until t'j_ reaches in its turn a crit- 
ical value where wj) signals a crossover toward a different 
instability of the metallic state where the Cooper pairing 
processes are prevailing and the instability is against s- 
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Fig. 4: Phonon frequency dependence of the superconducting 
ordering temperature Tc at different t'^ > t'^ and electron- 
phonon coupling gf^. The dashed lines correspond to the power 
law Tc ~ oJo"^. 



wave SC at Tc. According to Fig.[3]-a, the SC critical tem- 
perature increases with iojj and finally reaches a plateau 
above ujd ^ 2To, a scale apparently still tied to the adi- 
abatic CDW limit. The s-wave Cooper pairing attraction 
taking place along the chains and is strongly enhanced by 
CDW fluctuations (Fig.[l]-b). The T^ values thus achieved 
are markedly enhanced with respect to the BCS limit ob- 
tained when all the Peierls loops in ([^El) are put to zero. 

When in the adiabatic limit, t'^ is further increased up 
the critical value t'^ w 0.9To, nesting alterations are suf- 
ficiently large to suppress the singularity of the Peierls 
channel, bringing Tcdw down to zero (top curve. Fig. [S]- 
b). Since a,t lod — > 0, all the ladder diagrams are vanish- 
ingly small, the point t'^ defines a quantum critical point 
for CDW ordering. Moving now away from this point, 
along the frequency axis, the SC instability shows up as a 
result of the ladder diagrams of the Cooper channel that 
progressively unfold. However, at variance with the single 
channel BCS approximation, Tc follows a power law in- 
crease Tc ~ w^, with an exponent 77 ~ 0.70 smaller than 
unity (Fig.|4|) - the BCS value 77 ~ 1 being recovered when 
all the Peierls loops are turned off. 

Being independent of t'^ > t'j* and for a sizable range 
in weak coupling, the exponent 77 shows no noticeable 
trace of non universality (Fig. |4]). The non BCS increase 
of Tc is a direct consequence of the influence of CDW cor- 
relations on Cooper pairing and which takes place at all 
energy scales. As a matter of fact, when ojd increases, the 
Cooper diagrams grow in importance on one hand, but 
CDW correlations and then the pairing interaction is re- 
duced on the other. It is the combination of both effects 
that leads to an exponent 77 smaller than unity. This con- 
trasts with the BCS case where only the former effect is 
present, while the coupling is considered essentially fixed 
and attractive only below the sharp cutoff Eq{£)/2 = ujjj 
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The 'critical line' ^ lu^I) corresponds to a quantum- 
classical transition between the metallic and the SC states. 
The Debye frequency can thus be put in the category of 
a symmetry breaking parameter that drives the transition 
at a fixed clectron-phonon coupling. In the standard ter- 
minology [21| . ry is the crossover exponent (f> ~ ziy of the 
transition expressed in terms of the product of the dynam- 
ical and coherence length exponents. Assuming Lorcntz 
invariancc of the model, this forces z = 1, which would 
imply an anomalous dimension for the coherence length 
exponent, namely v = rj. Interestingly, if one looks at 
the range of phonon frequency over which the power law 
for Tc takes place, one realized that it is confined to low 
frequency. According to Figs. [2]-a and [21 Tc indeed levels 
off when ujd exceeds the Peierls scale 2Tq for non adia- 
baticity, stressing once again the non BCS character of 
the transition. 

To complete the analysis of the transition profile as a 
function of nesting deviations, one observes from Fig. [3]-b 
that at finite lod, the weakening of Tcdw by ladder dia- 
grams and vertex corrections is correlated to a reduction 
of the threshold value t'^ for the onset of superconductiv- 
ity. The decrease of Tc with t'^ > t'l is found in Fig. [3]-b 
to be relatively slow for any finite lud- In effect, t'j^ is more 
effective as a low energy scale to cut the CDW singularity 
than to CDW correlations responsible of the major part 
of the attractive pairing (Fig. ^h) . 

The above one-loop RG results for the present model 
are summed up in the global phase diagram of Figure [1] 
for arbitrary nesting deviations and phonon frequency. 

Connection with experiments. The above re- 
sults can apply to low dimensional CDW systems where 
the electron-phonon interaction is the prevalent mecha- 
nism for ordering. This is distinctly possible in CDW 
compounds like Per2[Au(mnt)2] for which the flat Pery- 
lene molecular unit is rather large in size and polarizable. 
The resulting Coulomb interaction is thus expected to be 
small as corroborated by the weak enhancement of the 
electron spin susceptibility and nuclear spin-lattice relax- 
ation rate reported for this material [22ll23] . The electron- 
phonon model considered above can then in a first approx- 
imation be applied to the CDW-SC sequence displayed by 
the compound under pressure. Taking into consideration 
the pronounced quasi-lD anisotropy of this material [3|, 
the TcDw(— 12K [IE]) observed in normal pressure condi- 
tions can be considered not too far below the optimal scale 
To calculated at perfect nesting. The Debye frequency for 
the 2kF acoustic phonons of the quarter-filled Perylenc 
stacks, though not known with accuracy, can be at least be 
taken as few dozens of degrees. This fairly places the com- 
pound with a frequency ratio uj]j/2Tq > 1, and according 
to Fig. [TJ with favorable conditions for superconductivity 
under pressure. As pressure scales up both band param- 
eters and phonon frequency [Mj, the system is likely to 
move from the CDW region {t'^ < t'l^uo /2T(^ > 1) to- 
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Fig. 5: Variation of the transition temperature (Tcdw: 
squares; Tc- triangles) as a function of the antinesting param- 
eter t'^ for different Coulomb interaction amplitude (7/|(j5''^|, 
normalized by [gf'l along U = 2V , and for ujd/2Ti) = 1. The 
dashed line gives the reduction of the maximum Tc with U . 
The circles (stars) refer to the SDW (d-wave SC) instability 
emerging for U /\g\^\ > 1. 



ward SC where [t'^ > I'I.ujdI'^Tq > 1) in the (^^,^0) 
plane of Fig. [TJ a path congruent with the results of Graf 
eA, al. [I]. 

The pairing attraction of the model, albeit boosted by 
CDW correlations, takes place along the chains and is re- 
sponsible for the s-wave character of superconductivity. 
While a s-wave order parameter is well known to sustain 
the presence of non magnetic impurities; on the contrary, 
it is sensitive to Coulomb interaction that is finite in prac- 
tice and acts as a pair breaking effect. The impact the 
Coulomb term has on the above results can be readily 
examined by modifying the initial conditions of the fiow 
equations and In the framework of the extended 
Hubbard model, this amounts to add the constant terms 
gi ~ U and g2 ^ U + 2V to the backward and forward 
scattering amplitudes at quarter-filling. Here U > and 
y > are the on-site and first nearest-neighbor repul- 
sive intra-chain interaction parameters of the extended- 
Hubbard model, here normalized by irvp. 

From the foregoing analysis and for suitable conditions 
for superconductivity at ujd/2To ~ 1, the one-loop RG so- 
lution for the critical temperature along the line U = 2V is 
given in Fig. [Sj as a function of the antinesting parameter 
t'^. While weak intra-chain Coulomb interaction reduces 
slightly Tcdw (essentially due to the repulsive backscat- 
tering component gi ) , its detrimental impact on supercon- 
ductivity is particularly pronounced for relatively weak re- 
pulsive interactions. The maximum Tc at drops by an 
order of magnitude at U ^ \gi^\/5. However, it worth 
noticing that such a range of U is sufficient to reduce 
the ratio Tc/TcDwImax to values comparable with the one 
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found experimentally in Per2[Au(nint)2] [I] and in the in- 
organic trichalcogenide compound NbSes [3- 

By reinforcing repulsive interactions, will therefore 
not be very long to become vanishingly small and po- 
tentially undetectable in practice, in spite of a sizable 
TcDW a-t low pressure. This may supply some insight as 
to why quarter-filled molecular compounds like TMTSF- 
DMTCNQ gg, known as a correlated CDW system pi^ . 
failed to show any sign of superconductivity following 
the suppression of its lattice distorted state under pres- 
sure. Such an interpretation may also adhere to the ab- 
sence of superconductivity in the phase diagram of cor- 
related quasi-lD compounds like [EDT-TTF-C0NMe2]X 
[17], and (DI-DCNQljsX [H] at high pressure. From this 
angle, the chance for a correlated organic metal like TTF- 
TCNQ [S] , whose CDW order is expected to vanish around 
90 kbar [29], to show the presence of superconductivity is 
reduced at the very least in the s-wave channel. 

We finally examine in Fig. [S] the impact of further 
increasing the Coulomb term U on the phase diagram, 
namely beyond the amplitude of the phonon induced in- 
teraction. A qualitative change then occurs in the nature 
and the sequence of ground states. The CDW gives way to 
a SDW instability at low t^, which in its turn yields non 
conventional d-wave superconductivity at t'^ and above. 
It is where the results of the present work connect to the 
sequence of instabilities found in the purely repulsive case 
[Tmn] , which is known to apply in systems like the Bech- 
gaard salts [5HTU]. 

Conclusion. — In conclusion the instability of the 
metallic state against charge-density-wave and s-wave su- 
perconducting orders in quasi-lD systems can be analyzed 
by extending the RG approach to electrons coupled to 
phonons of arbitrary frequency. The results show that 
both instabilities influence one another and form a se- 
quence of ordered states that captures the key traits of 
low dimensional charge density- wave materials exhibiting 
superconductivity under pressure. 
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